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The structural simplicity of two-terminal compared to three-terminal 
devices indicates the potential importance of two terminal devices employing 
semiconductors and having negative resistance at frequencies properly re- 
lated to the transit time of carriers through them. Such negative resistances 
may be combined with un symmetrically transmitting components, such as 
gyrators or Hall effect plates, to form dissected amplifiers that may be made 
to simulate conventional three-terminal amplifiers and operate at high fre- 
quencies. The characteristics of several structures are analyzed on the basis 
of theory and it is found that negative resistances are possible for properly 
designed structures. 

1. NEGATIVE RESISTANCE AND DISSECTED AMPLIFIERS 

Because the drift velocities of current carriers in semiconductors are 
smaller than the velocities attainable in vacuum tubes, transistor struc- 
tures must be smaller to achieve comparable frequencies. In principle 
it is possible, of course, to make compositional structures (i.e., distribu- 
tions of donors and acceptors) in semiconductor crystals on a scale much 
smaller than is possible for vacuum tubes. At present, however, the 
available techniques are limited and it may require many years before 
the ultimate potentialities are approached. 
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It is instructive, however, to speculate on some of these ultimate po- 
tentialities. For example a grain boundary formed of edge type disloca- 
tions is in a sense an analogue of a grid. Possibly it can be made into a 
grid by acting as a locus for an atmosphere of donors or acceptors. 
Evidently such a grid will approach the smallest spacing that can be 
achieved with any known form of matter. If the spacings perpendicular 
to the grid are made comparable to a mean free path of the carriers used, 
the device will operate like a vacuum tube with carrier velocities con- 
trolled by inertia rather than by mobility. It is not easy to conceive of a 
structure having the potentiality of operating at higher frequencies. 

It is evident that the difficulty of making small structures increases 
with the number of electrodes. For example, it is now possible to make 
diodes which give usable rectification at frequencies above 10 10 cps. In 
these the "working volume" is a very thin layer under the metal point. 
The thickness of this layer is controlled by surface treatments and the 
applied voltages. The diameter of the point, which is the mi n imum di- 
mension mechanically controlled, is much larger than this thickness of 
the layer. In order to make a transistor of comparable frequency, it 
would be necessary to make structural elements having dimensions 
comparable to the thickness of the layer and this would be a much more 
exacting task than making the diode. 

These considerations point out the importance of giving serious con- 
sideration to two-terminal structures as amplifying elements. It is pos- 
sible, in principle at least, to have structures which are much smaller 
in one dimension than the other two and which exhibit negative re- 
sistance and thus give ac power at frequencies comparable to the 
reciprocal of the transit time across the small dimension. 

The attractiveness of such negative resistance diodes for amplification 
is enhanced by the possibility of using them in dissected amplifiers ' in 
combination with nonreciprocal elements such as gyrators or Hall effect 
plates. Combinations of negative resistance elements and nonreciprocal 
elements can lead to structures having gain and unsymmetrical trans- 
mission that simulate conventional amplifiers. The adjective dissected 
has been suggested for them since elements giving power gain are 
physically separated from those giving one-way transmission. 

In this article we shall not consider the possible forms of dissected 
amplifiers, of which there are a wide variety. Instead we shall give an 
introductory treatment of some forms of negative resistance that may 
arise from transit time effects. In some cases the most instructive way 
of treating the structure is by way of the "impulsive impedance" and 
we devote most of the next section to considering this method. 
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2. THE IMPULSIVE IMPEDANCE AND NEGATIVE RESISTANCE 

The impulsive impedance D(t) for a two terminal device is denned in 
terms of its transient response to an impulse of current. Thus if the 
current through the device is 

/(*) = J + M, 
where J is the dc current and 

m = o 

except very near t = and 

f j(t)dt = 8Q, 

then the voltage is 

V(t) = V + v(t), 

where V is the dc voltage and 

v(t) = 8Q D(t). 

In other words, if in addition to the dc biasing current, a charge 8Q is 
instantaneously forced through the circuit at time t = 0, the added 
voltage is D(t) per unit charge. These equations also serve to introduce 
the notation used in this article: 

Notation. In general, quantities that are functions of time or position 
will have the functional dependence explicitly indicated. In Sections 4 
and 5, however, the symbol 5 will be used to distinguish the transient 
parts 8E and 8p from the dc parts of the electric field and charge density. 

Capital V(t) and J(t) stand for total voltages and current. Without 
functional dependence upon (t) they are the dc parts. Similarly v(t) and 
j(t) are the ac or transient parts. A sinusoidal disturbance is represented 
by 

v(t) = v exp iut, 

j(t) = J exp iwt. 

where v and j are not functions of time. Where it is necessary to distin- 
guish the displacement current at a particular location from the con- 
duction current, as in the next section, we shall write 
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meaning the displacement current across space charge region S 2 as a 
function of time. 

In this section we shall treat J andj as circuit currents. In subsequent 
sections, we shall be concerned with current densities and shall use the 
same symbols. 

The complex impedance of the device is evidently 

Z(a) = v/j, 

where v and j are the coefficients in the sinusoidal case. 

In terms of the system of notation introduced above, Z(u>) may also 
be expressed in terms of D(t) by expressing j exp iut in terms, of incre- 
ments of charge 

dQ = je*"* dt, 

and summing over all increments up to time t. This leads to 



Z(«) = [ D(t) exp {-wt)dt. 
Jo 

A negative resistance will occur if 

»oo *oo 

> / D(t) cos tatdt = (-1/co) / D\t) sin at dt, 

Jo Jo 

the latter form coming from integration by parts for the case of Z>(*>) 
= 0, the only situation treated in this article. 

The use of D(t) in analysing the potential merits of diode structures 
from the point of view of negative resistance is illustrated in Fig. 2.1. 
Here three cases of D(t) together with certain cosine waves are shown. 
It is seen for case (a) that a negative real part of Z will be obtained. 
For case (b), the real part of Z is zero for the frequency shown; this 
represents a limit; for other frequencies, a positive real part will be ob- 
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Fig. 2.1 — Some hypothetical D(t) characteristics. 
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tained. Case (c) represents an exponential fall such as might occur for a 
capacitor and resistor in parallel. We shall discuss this example below. 

The conclusions regarding (a) and (b) may be somewhat more easily 
seen from the corresponding —D'(t) plots shown in Fig. 2.2. From part 
(a) it can be seen that the negative maximum in the sine wave at the 
end of the rectangular D(t) plot is particularly favorable. From part (b) 
it is seen that no choice of o> will result in more negative area of sine 
wave than positive. For (c) it is evident that each positive half cycle of 
the sine wave gives a larger contribution than the following negative 
half cycle and hence that a positive resistance will be obtained. 

For case (c), it is instructive to obtain the value of Z analytically by 
using 



D(t) = C 1 exp (-t/RC). 



This leads correctly to 



Z(«) = (/T 1 + io>CT\ 



For small values of uRC, Z reduces to 7?; furthermore, for this case, the 
decay of D(t) occurs while cos at = 1. Under these conditions 



ZM = f 

Jo 



D(t) dt. 



This result is useful for estimating the effect of quickly decaying contri- 
butions to D(t). These evidently contribute a positive resistance to Z 
equal to the area under the D(t) curve. 

From these considerations it follows that an upward deviation from 
the linear fall in Fig. 2.1(b) towards Fig. 2.1(a) will result in negative 
resistance. In Sections 4 and 5 we shall see how particular structures 
may lead to such favorable, con vex-upwards characteristics for D(t). 
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Fig. 2.2 — The —D'(t) characteristic corresponding to fig. 2.1. 
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3. MINORITY CARRIER DELAY DIODE 

As a first example we shall consider the behavior of the device shown 
in Fig. 3.1. We have chosen a p-n-p structure rather than an n-p-n so as 
to deal with positively charged carriers and thus avoid numerous minus 
signs in the equations. In this figure we have used capital letters P and N 
to designate specific regions, reserving the small letters to indicate carrier 
densities and conductivity types. 

Several features that simplify the theoretical treatment should be 
noted : 

(a) The PiN junction is 100-fold more heavily doped on the Pj-side. 

(b) The doping in the layer N varies exponentially with distance by 
a factor of 10 across the layer. 
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Fig. 3.1 — Constitution of minority carrier delay diode. 






NEGATIVE RESISTANCE IN SEMICONDUCTOR DIODES 805 

(c) Throughout N the concentration of holes is less by a factor of 10 
than the electron concentration. 

(d) The thickness of N is large compared to the depth of space charge 
penetration into it. 

(e) The voltage drop across the space charge region S 2 is large com- 
pared to the other voltage drops. 

The conditions lead at the operating frequency to the following conse- 
quences : 

(1) The current across the first junction is carried preponderantly by 
holes. 

(2) The hole drift in N is substantially unaffected by the ac field and 
thus represents the delayed diffusing and drifting current injected across 
the first junction. 

(3) The ac voltage drop occurs chiefly across S 2 . 

We shall show below (Section 3.2) how (a) to (e) lead to (1) to (3), 
but we shall first bring out the importance of (1) to (3) by using them to 
give a simple treatment of the impedance of the diode. 

3.1. Calculation of Impedance. 

If the total current is 

J(t) = J+je iut , (3.1) 

then the ac hole current across Si is also in the notation discussed in 
Section 2 with the addition of the symbol p to indicate holes 

j(p, Si, t) = je iut . (3.2) 

This current flows through the n-layer unaffected by the ac field and 
arrives at <S 2 delayed and attenuated by a complex factor 

/3 = | | exp (-10). (3.3) 

Because of the high field in *S> 2 , the transit time there is negligible so that 
the hole current arriving at P 2 is 

j(p, Sfr M = &<"*. (3.4) 

In addition to this current, there is a dielectric displacement current in 
S 2 which is converted to hole conduction current in P 2 . If the voltage 
drop across S 2 is 

V(S 2 , t) = 7(S,) + v(S 2 )e iut , (3.5) 
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then the ac displacement current is 

j(D, S 2 )e iut = uaC#(S 1 )e i " t . (3.6) 

Now the total current is constant through the device, hence 

3 = 3(V, St) + J(D, &) (3.7) 

which leads to 

j = faC*(A)/(l - &). (3.8) 

If v(S 2 ) is substantially equal to the ac voltage across the unit, then 
the impedance is 

Z = v(Si)/j = iX/iuCi) + ifi/aCt 

= (1A«C 2 ) + (l/wC 2 ) 1 18 | exp t[(ir/2) - 09. 

Evidently if > t and 6 < 2x, the second term will have a negative 
real part so that the diode will act as a power source. 

If we neglect the ac electric field in ./V then § may be calculated in 
terms of the thickness L = x 2 — xi of the layer and the potential drop 
across the layer. This latter arises from the concentration ratio N d \/N d2 
between the two sides of N. Since the donor charge density is neutralized 
substantially entirely by electrons, and since almost no electron current 
flows, the electron concentration difference must result from a Boltzman 
factor (at 10 17 /cm 3 Fermi-Dirac statistics are not needed) and this leads 
to 

AV n = (kT/q)ln(N dl /N d 2) (3.10) 

for the potential drop across N. In N the electric field is thus 

E = AF„/L. (3.11) 

The differential equation for hole concentration for a disturbance of 
frequency a> is 

J >-i B p-_ M »B*2 + D,g. (3.12) 

This linear differential equation has two linearly independent solutions. 
These must satisfy at x 2 , the left edge of the space charge layer £ 2 , the 
boundary condition that the hole density is practically zero. The 
appropriate solution is 



■p = e 



iut r Ai(i-x 2 ) _*2<*—*2> 



- e k * { *- X2) ]. (3.13) 
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where 



where 



Lh m (xi - xi)h = «[1 + (1 + iy) m ), 
Lk 2 m (x 2 - Xi)h = a[l - (1 + iy) m ] 



a = qAV/2kT, 

y = 4wD p /w 2 , 

u = n P E = ti p AV/L. 



The current is 



(3.14) 
(3.15) 

(3.16) 
(3.17) 
(3.18) 

(3.19) 



j(p, x, t) = q(up - D p dp/dx) 

and the ratio of currents at X\ and x% , which is /3 by definition, is 

(3 = j(p, x-i , l)/(j(p, xi , I), 

Lik\ — L/ko 
~ Kh exp (-Lk 2 ) - Kk 2 exp (-LkJ (3.20) 

= 2(1 + iy) m e a 

[1 + (1 +z'7) ,/2 ]expa(l +*y) 1 « - [1 - (l+n) 1/8 ]exp - ail+iy) 1 ' 2 ' 

The phase lag in /3 mvist exceed 180° or it in order to give negative re- 
sistance. It can be seen that this phase factor must result from the first 
exponential in the denominator by the line of reasoning suggested below: 
The real part of the exponent is larger than the imaginary part. Hence 
the absolute ratio of the two exponentials is at least 2x. For this condi- 
tion the second term in the denominator is negligible compared to the 
first. Hence the phase of (3.20) is determined largely by the first ex- 
ponential. As a helpful approximation we may neglect the second term 
and write 



2(1 + iy) m exp [a - a (1 + iy) m ] 

P 1 + (l + ty) 112 

Two limiting cases are worthy of special note: 

(I) a — > 0, uniform n-layer, y — » °o . 

= 2 exp -a(iy) m = 2 exp -(1 + i){w/2D) m L. 

(II) a-> oo,qAV/kTy> l,T-*0. 
= exp[-*ay/2 - ay 2 /8], 

= exp[-i(uL/u) - (DL/w)/(m/ w ) 2 J. 



(3.21) 



(3.22) 



(3.23) 
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These expressions may be interpreted as follows: In case (I), flow is by 
diffusion and the propagation factors fa and k% take the form 
±(a/L)(iy) m . For this case attenuation and delay terms in the ex- 
ponential are equal, and the largest negative term occurs in Z when the 
phase angle is 5tt/4 (as may be verified by differentiation.) This leads to 

- 2(-l + i)2~ m exp (-5tt/4) = 0.028(-l + i), (3.24) 

which gives 

Z = (l/oAX- 0.028 -i 1.028), (3.25) 

the impedance of a condenser with a negative Q of 37. In order to make 
an oscillator by coupling this to an inductance, an inductance with a Q 
of more than 37 must be used. It is obviously advantageous to reduce 
the magnitude of the negative Q. 

For case (II) in its ideal form, the ac current simply drifts through the 
n-layer without attenuation. This produces a phase lag of u times the 
transit time L/u. If this were the only effect involved, a capacitor with 
a negative Q of less than unity could be produced. In addition, however, 
there is attenuation due to spreading by diffusion. This effect is depend- 
ent upon the ratio of the spread by diffusion (DL/u) m to the separation 
of planes of equal phase in the drifting hole current. This latter separa- 
tion is 2iru/io. The square of this ratio appears in the attenuation term 
in the second form of /3. 

We shall estimate the effect of the attenuation term by taking 

cry/2 = 3tt/2, (3.26) 

so that the desired phase shift is obtained. The attenuation term is 
7/4 smaller than this so that if t/4 is considerably less than one, the 
attenuation in will be small while the phase shift is correct. If we take 
3t/2 for the value of ay/2, then the value of 7 becomes 

7 = 4o>D/u 2 = QTkT/qAV. (3.27) 

Thus the approximation on which (II) is based fails unless q&V/kT 
> 18, a value which implies an enormous range of concentration in the 
n-layer. We must, therefore, investigate the case of gradients in the 
n-layer by more complete algebraic procedures. 

We shall denote by — 0i the phase shift in due to the exponential in 
equation (3.21). The total phase shift is somewhat less since the alge- 
braic expressions give a small positive phase shift of at most about 15°, 
which vanishes for large and small values of 7. Similarly the attenuation 
of arises chiefly from the real part of the exponential since the absolute 






NEGATIVE RESISTANCE IN SEMICONDUCTOR DIODES 809 

value of the algebraic expressions lies between 2 for 7 = and 1 for 

y = 00 . 

It is instructive to express the real part of the exponent in terms of 
a and 6\ . This is done as follows: 

a - a(l + iy) m = -17 - iBx . (3.28) 

This can be solved for 77 and (1 + iy) m : 

r, = (a 2 + 0f) 1/2 - a, (3.29) 

(1 + iy) m = [(a 2 + 0! 2 ) 1 ' 2 + id 1 ]/a. (3.30) 

From there it is seen that for a fixed value of 0i , the attenuation can be 
greatly reduced by increasing a. Unfortunately, this requires very large 
changes in concentration. For example with 61 = 3t/2 and a = 81 , the 
value of 77 is reduced to 77 = 0.414 0i . However, the value of potential 
difference is 

qAV/kT = 3tt, (3.31) 

giving 

N dl /N d2 = 10\ (3.32) 

For the case shown in Fig. 3.1, a = 1.15 and for 6\ = 57r/4 = 3.9 we 
obtain 

77 = 2.95 = ln e 19.5 (3.33) 

This is an improvement of about 1 factor of e in the exponential com- 
pared to having AV = 0. The value of /3 is 

= 0.082 X exp (-i 218°), (3.34) 

and this leads to 

Z = (coC 2 ) _1 (-0.05 - % 1.065). (3.35) 

Thus at the operating frequency, the diode appears to be a capacitor 
with a negative Q of 21. 

Increasing the concentration change to a factor of 100, so that a = 
2.3, gives 

77 = 2.25, (3.36) 

/3 = 0.18 < - 220°, (3.37) 

Z = (aA)" 1 (-0.116 - i 1. 14), (3.38) 
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Q = -10. (3.39) 

The calculations indicate that attenuation can be controlled to a 
considerable degree while maintaining the desired phase shift. 

3.2. Justification of Consequences (1), (2) and (3) 

In germanium at room temperature the product np is about 10 
under equilibrium conditions. At the first junction of Fig. 3.1 it is 10 , 
implying a forward bias 3 ' 2 of (kT/q) 2.3 X 5. In order to maintain this 
forward bias a flow of electrons must be furnished to N. There are sev- 
eral ways of accomplishing this. In the first place, the reverse bias 
across & draws a reverse current of thermally generated electrons from 
jP 2 . This current can be controlled by controlling the lifetime and 
temperature in the P 2 region. Alternatively, electrons may be injected 
into Pi ; some of these will diffuse to £ 2 and arrive at N. Still another 
means of controlling the bias across & is to make contact to N itself. 
Since only the dc bias need be controlled, the series resistance across N 
itself is unimportant; the source should be of high impedance. 

The decrease in density of 10 4 across the junction in carrier concentra- 
tion implies a potential difference of 9.2 (kT/q). Most of this potential 
difference occurs where the carrier concentration is negligible. Hence 
the space charge theory may be applied. Furthermore, the acceptor con- 
centration is much higher than the donor concentration. Hence the 
space charge extends chiefly into the donor region and we may write ' 

AFi = (2irqN d /K)W 2 (3.40) 

for the relationship between width W of the space charge region and 
voltage drop AV. 

If this voltage drop has an ac component, then a charging current will 
be required to change W. This current is determined by the admittance 

uC = uk/4htW (3.41) 

of the space charge region. 

At the same time injected hole and electron currents flow across the 
junction. The admittance associated with the hole current is approxi- 
mately 

A = (ia/D) m a pl , (3.42) 

where <x pX is the hole conductivity just inside the n-layer. 3 ' 4 Actually, as 
discussed below, the admittance is somewhat higher. 




NEGATIVE RESISTANCE IN SEMICONDUCTOR DIODES 811 

The ratio of the admittances is 



(3.43) 



For our example this expression is much greater than 1 as may be seen 
as follows : The first fraction is the ratio of the dielectric decay constant 
to to. This is 10 3 or more larger than a> need be. The next term is about 10 
and the last term is the ratio of hole to electron density at Xi and is 
about 10 -2 . Hence the ratio of impedances is about 10: 1. 

We shall next consider why the expression for A for holes must be 
examined more closely. The admittance formula used above applies to 
the case of zero field to the right of the junction. The aiding field will 
increase the flow of holes into the n-layer and raise the admittance 
somewhat. Correcting for this will increase A in respect to coC and will 
thus strengthen rather than weaken the argument. 

Also in the expression for A, no account was taken of the transit time 
across the region W. If we assume a uniform field in this region for 
purposes of making estimates, then the solution of equation (3.20) may 
be applied. Since now u corresponds to drift velocity due to QkT/q of 
voltage drop across W which is much less than L in length, it is evident 
that 7 will be less by a large factor in this region compared to its value 
in N. This leads to the conclusion that phase lags will be unimportant 
in this region. 

We have neglected the effect of electron injection into Pi . By the 
customary arguments for unsymmetrically doped junctions, it follows 
that this current is very small compared to the hole current. 

This justifies consequence (1). 

Consequence (2) may be justified as follows: At Xi all the ac current 
j exp (iut) is carried by holes. If a pure drift case occurred, the hole cur- 
rent might be reversed at some point in the n-layer and be —j exp (iut). 
Under these conditions the electron current would have to be 2j exp 
(iut). Under no conditions, however, will the electron current be larger 
than this. This maximum possible electron current will require an electric 
field and this field will also affect the hole flow. Since the electron con- 
ductivity is at least 10 times larger than the hole conductivity, the hole 
current due to the ac field will only be about ^ o of j at most. Thus the 
hole current is only slightly affected by the ac field. 

Consequence (3) follows from the fact that the reverse biased junction 
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£2 has much higher impedance than Si . & has higher impedance than 
that for hole injection into N. The impedance for hole injection into N 
corresponds to the hole conductivity in the n-layer over a distance com- 
parable with the thickness of N. However, the impedance of N itself is 
that due to the much larger number of electrons in it and is thus much 
less than the impedance of 81 . Thus it follows that the impedances 
across & and across N are much less than across £ 2 • This conclusion is 
not affected by the modification of impedance of S 2 due to hole flow 
across it. 

3.3. Modifications 

The treatment presented above has been based upon the conditions 
(a) to (e). Some of these are advantageous from the point of view of 
operation but others have been introduced to simplify the treatment. 
Among the latter is the condition that the current across 81 is carried 
chiefly by holes. If the current were chiefly capacitative at this junction, 
then the voltage would lag 90° behind the current. This adds a desirable 
phase lag in the hole injection across Si and thus requires less phase shift 
in the n-layer. By suitably adjusting the ratio of capacitative and in- 
ductive admittances, a net improvement in Q may be obtained. 

4. THE TRANSIENT RESPONSE IN A UNIPOLAR STRUCTURE 

In the previous section the electric field produced by the injected holes 
had a negligible influence on the motion of the injected holes. In effect 



(a) 




LU 0. 



(b) 

Fig. 4.1 — Spnce charge limited hole flow. 
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this was due to the bipolar nature of the mode of operation considered, 
the majority carriers in the region N acting to shield the minority car- 
riers from their own space charge. 

In this section we shall deal with unipolar diodes in which only one 
type of carrier is present in sufficient number to have a major effect. In 
these the influence of the space charge of the carriers upon their motion 
plays an important role. 

Fig. 4.1 illustrates one example of the type of structure covered by 
the theory of this section. It is again a p-n-p structure like that considered 
in Section 3. However, in this case the dimensions, the donor density 
and the applied potential are such that the space charge "punches 
through" the device. 4 ' 1 Under these circumstances a condition of space 
charge limited emission is set up so that holes are injected from the 
positive region Pi to just such an extent that their flow is limited by 
their own space charge. This limitation is associated with the maximum 
of potential just inside N. 

The potential maximum is evidently a "hook" for electrons generated 
thermally in P 2 and in N. Under some circumstances electrons may 
accumulate and form a layer in which there is no electron flow and hole 
flow is carried equally by diffusion and drift. Such stagnant regions will 
tend to be suppressed if Pi is made of short lifetime material, so that 
electrons are siphoned out of N, or if p at the maximum is larger than p 
for intrinsic material and the lifetime is locally low. 

We shall treat the transient response of this structure of Fig. 4.1 from 
the point of view of the impulsive impedance discussed in Section 2. 
Accordingly we suppose that a steady current J flows per unit area. At 
t = an added pulse of current occurs carrying a total charge of 8Qi per 
unit area, the subscript "i" signifying initial condition. Our problem is 
to determine how this added charge is carried by a transient disturbance 
in the hole flow and what is the resultant dependence of voltage upon 
time ; by definition the added voltage across the device is 

v(t) = 8QiD(t). (4.1) 

Since we are dealing with a planar model, we shall suppose that the 
initial condition at t = corresponds to added charges 8Qi and — 8Qi on 
the metal plates on the P-regions. These charges set up an added field 

8E { = 8Qi/K, (4.2) 

where 

K = k€ (4.3) 



. 
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in MKS units. The initial value t>(0) is then simply 8E { times the total 
width of the structure. 

The first effect, which takes place in a negligible time in respect to 
the frequencies involved, is the dielectric relaxation of the field in Pi 
and P 2 . The added current due to 5/?, leads to an exponential decay of 
8E in these regions with a transfer of 5Q> and — 8Q> to the two boundaries 
of iV. If Pi and P 2 are thin compared to JV, the resulting drop in v(t) is 
small. In any event it can be shown by the reasoning at the end of Sec- 
tion 2 that this contribution to D(l) adds simply the series resistance of 
Pi and Pi to the impedance. 

The next effect is the transport of 5Q, on left side into N by hole flow 
over the potential maximum. It will be easier, however, to discuss this 
process after the treatment of the transient effects that occur in N 
itself. Consequently, we shall at this point assume that after a time, 
short compared with the important relaxation time in the structure, the 
disturbance of hole density is as shown in Fig. 4.2(a). 

Fig. 4.2(a) shows added charges +5Q; and -5Q, produced by a 
disturbance denoted as 8p in the hole density. The charge -5Q, on the 
right side is produced by an increased penetration of the space charge 
into P 2 ; it is similar to that produced by increasing reverse bias on a p-n 
junction. 

Fig. 4.2(b) shows the corresponding disturbance in electric field. This 
disturbance is denoted by 8E which is a function of x and t. Evidently 



v{t) = f 
Jo 



8E(x, t) dx. 



(4.4) 



and this is the area under the 8E curve. 

The other parts of the figure indicate qualitatively a subsequent stage 
in the motion and decay of 8p and 8E. Our problem is to formulate mathe- 
matically this decay process. We shall treat the decay process in terms 
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Fig. 4.2 — The initial stage and a subsequent stage of the transient. 
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of the effect of drift -in the electric field and neglect the effects of diffu- 
sion. This procedure can be justified by the fact that as soon as a hole 
had reached a point where the potential has fallen by kT/q below the 
maximum, its flow is governed by drift rather than diffusion and the 
predominance of drift continues to increase towards the right. 4 ' 2 

If drift in the field is the predominant cause of hole flow, then the 
equations governing the situation in N are 

J + 8 J = (p + 8p)(u + 8u), (4.5) 

where the terms with 5 represent the transient effects and those without 
represent the steady state solution, p = qp is the charge density of the 
holes and u their drift velocity. The equation for the change of E with 
distance is 

K(d/dx)(E + 8E) = Pf + p + 8 P , (4.6) 

where p/ is the fixed charge density due to donors and acceptors. (We 
neglect any effect of traps.) The steady state equation for E is thus 

K(dE/dx) = p, + J/u. (4.7) 

In a region where p f is independent of x, this equation may be reduced 
to quadratures by writing 

K dE/( P{ + J/u) = dx; (4.8) 

the left side is then a known function of E through the dependence of 
u upon E. 

It is convenient to introduce a time-like variable s which is the transit 
time for the dc solution. Evidently 

ds = dx/u = KdE/(p f u + J). (4.9) 

For the case of space charge limited current, s may be conveniently 
measured from the potential maximum. Even though the solution is 
invalid at that point, the integrals converge and the contribution from 
the region within kT/q of the maximum is small. 

We shall assume that the equations for the steady state case have been 
solved and that the functional relationships are known between E, x, v 
and s. 

The differential equation for 8E may then be obtained as follows: 
To the left of the pulse in 8p in Fig. 4.2(a), 8E is zero. From equation 
(4.6) we have 

KdE(x)/dx - 8 P . (4.10) 



. 









(4.11) 
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Integrating this from the region where E is zero gives 

K8E(x, t) = [ 8 P (x, t) dx. 
Jo 

Equation (4.11) states that the dielectric displacement at x is equal to 
the excess charge between the potential maximum and x. Evidently 
during the transient following Fig. 4.2(a), the rate of change of this 
extra charge is — 8J(x, t) since the dc current is flowing in at the left 
and an excess current 8J flows out at the right. Hence we have 






(4.12) 



Kd8E/dt = -8J, 

= — (8pu + p8u). 

For the change in drift velocity we may write 

8u = (du/dE) 8E = p*8E. (4.13) 

For high electric fields u increases less rapidly than linearly with E and 
p* is less than the low-field mobility. 4 3 For very high fields p* is nearly 
zero and there are theoretical reasons for thinking that there may be a 
range in which p* is negative. We shall return to this point in the next 
section. 

In Fig. 4.3 we show a diagrammatic representation of the transient so- 
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Fig. 4.3 — Graphical representation of the dependence of dE upon time. 
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lution. Each of the dashed lines represents the decay of 8E as measured in 
a moving coordinate system : Thus we consider 8E measured at a position 
x(«a + 0; this is a position that moves with the dc velocity u. This 8E 
is evidently expressed in terms of 8E(x, t) by writing x = x(so + I): 

8E in moving system ■ 8E m (s , t) = 8E[x(s + t),t]. (4.14) 

The differential equation for 8E m is 

(d/dt) 8E m = (d8E/dt) x + (d8E/dx) t dx/dt, 

= (d8E/dt) x + (d8E/dx) t u, 

(4.15) 
= -(u8 P + P 8u)/K + {8p/K)u, 

= -(pn*/K)8E = -v8E, 

where the quantity 

v = pp*/m (4.16) 

is an effective dielectric relaxation constant being the differential con- 
ductivity pp.* divided by the permittivity K. 

Evidently v is a function of position x only and may be expressed as 
v(s) through the dependence of x upon s. Thus we may write 

(d/dt)8E m (so , t) = -p(s + t)8E m (so , t) (4.17) 

which has a solution 

8E n (s , t) = 8E m (s , 0) exp [-g(s + t) + g(s )], (4.18) 

where 

p«o+t 

g(s + t) = v(s) ds. (4.19) 

The lower limit s' is chosen for convenience so as to avoid singularities in 
g(s). This integration shows that 8E m decays exponentially as the elec- 
trical field would decay in a material whose dielectric relaxation constant 
changed with time just as v changes as observed on the moving plane. 

Fig. 4.3 shows on the dashed lines the decay of 8E m on the moving 
planes. Since 8E m is zero to the left of the initial pulse in Fig. 4.2(a), it 
remains zero on all moving planes which follow the pulse of 8Qo . This 
justifies the statement made earlier. The solid curves labelled ti , h etc. 
show the spatial dependence of 8E for times h , U , etc. after the charge 
8Qi is added. 

The values of the transient voltage v(t) at time h , for example, is the 
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integral under the curve k . This curve is zero for x < x(h) and for 
x > x(ti) it is 

8E(x, h) = (SQi/K) exp [- g(s + U) + g(s )], (4.20) 

where 

x = x(s + 0- ( 4 - 21 ) 

If the total transit time across N is S so that 

x(S) = L, (4.22) 

then 

v(h) = [ 8E(x, h) dx. (4.23) 

From this expression we can derive the desired formula for D(t). For 
this purpose the integral over dx is replaced by an integral over s. At 
time t the range of s is evidently from t to S and dx = u(s) ds. From 
this we obtain: 

D(t) = v(t)/5Q i} 

-i (4.24) 

= (1/10 ] exp [- g(s) + g(s - t)]v(s) ds. 

From Fig. 4.3 we can see that there are competing tendencies in the 
decay of D(t) some of which tend to produce the desired convex shape 
discussed in Section 2 and others the concave shape. The effect of the 
dielectric relaxation constant is adverse and tends to produce an ex- 
ponential decay. On the other hand the advance of the pulse of holes 
from left to right in Fig. 4.2 proceeds in an accelerated fashion with the 
result that the range of x over which 8E is not zero decreases at an ac- 
celerated rate. If the dielectric relaxation were zero, this would result 
in the desired convex upwards shape. 

The resultant shape of the D(t) curve is thus sensitive to the exact 
relationship of the transit time and dielectric relaxation. This can be 
illustrated by giving the results of analysis for a p-n-p structure, neglect- 
ing diffusion and considering n to be constant. The solutions of the dc 
equations are readily obtained for this case and have been published. 4 4 
For convenience we repeat them here: 

E = (J/w)(e a ° - 1), (4.25) 

x(s) = nJK{np f T 2 (e al -at- 1), (4.26) 









NEGATIVE RESISTANCE IN SEMICONDUCTOR DIODES 



819 



L = x(s) = (JK/ m 2 )(<f - - 1), 



= as. 



From these it is found that 



This leads to 



lng(s) = (1 - e~ aa ). 
D(t) = (J /up?) \i + e al (at- 0- 1)] 

EE (//MP/) D(/3, O0. 



(4.27) 



(4.29) 



(4.30) 



For * = this reduces correctly to L/K. 

Figure 4.4 shows the resulting shape of the D curves with /3 as a 
parameter. Large values of /3 correspond to cases in which the hole 
charge density is small compared to p/ and to relatively long relaxation 
constants. For them the desired convex upward shape results. 

Figure 4.5(a) and 4.5(b) show the real and imaginary parts of the 
impedance expressed in terms of Z(/3, 6) : 



Z(») = [ ' 

Jo 



e~ iai D{t) dt 



= (Kj/nV f )z(p, e\ 

= u>T. 



(4.31) 



(4.32) 



O 8 

o 

s 

Q 0.6 
\ 

s 

Q 0.4 



0.2 

























8 s V\ 


>£ X 


^ V 

























"0 0.2 0.4 0.6 0.8 1.0 

cti//3 
Fig. 4.4 — Impulsive impedance for various values of in p-n-p structure. 
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It is seen that values of — Q as small as about 10 can be obtained for 

^ 3. 

In the next section we shall consider modifications which may result 
from variations in n* and for changes in geometry. 

We must return to the question of how the charge +8Q t passes the 
potential maximum. In order that the theory given above apply, it is 
necessary that the time required for 8Qi to enter the drift region be short 
compared to the transit time. At the potential maximum the charge 
density may be estimated by the methods previously dealt with in the 
theory of space charge limited emission. Initially +SQi appears to the 
left of the maximum and the field at the maximum is 8E { . This field will 
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Fig. 4.5 — Impedance of a p-n-p structure, (c) Real part of impedance. 
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then relax with a relaxation constant of about fip(max)/K where p(max) 
is the hole charge density. Actually the relaxation may be somewhat 
quicker because the concentration gradient of the added holes also 
contributes to the flow over the maximum. Since the charge density at 
kT/q below the maximum is comparable to that at the maximum the 
entire relaxation process will proceed at about this rate. Thus a criterion 
for the applicability of the theory is that K/np(max) be much less than 
S, the transit time or total decay time for D(t). 

5. MOBILITY AND GEOMETRY EFFECTS 

5.1. The Effect of a Region of Negative /** 

In very high electric fields holes may be expected on the basis of 
theory to exhibit a negative value of n*. This theory 6 ' 1 is founded on the 
idea that a hole can lose energy to phonons at a certain maximum aver- 
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age rate P ma x . (Pmax is the staukonstante of Kromer.) The holes probably 
achieve this rate when their energy is near the middle of the valence 
band. Under these conditions the power input from the electric field 
must be no greater than P max : 



From this it follows that 



qEu g P, 



u ^ P m **/qE, 



(5.1) 



(5.2) 



so that the drift velocity will decrease with increasing field at sufficiently 
high fields. 

Furthermore, if the width of the valence band is less than the energy 
gap, then a hole cannot acquire enough energy to produce hole electron 
pairs. Thus in such a case, the negative resistance range should be reached 
before breakdown effects occur. 

In Fig. 5.1 we illustrate the general trends of the u versus E curve, to 
be expected if the staueffekt occurs. As is indicated, the maximum drift 
velocity will be referred to as u m . It occurs at a field E m . Since we are 
here concerned with principles rather than details, no attempt has been 
made to indicate the square root range in which u is proportional to E . 
This range has been observed by E. J. Ryder 6,2 and shown by G. C. 
Dacey 5 " 3 to control hole flow in space charge limited hole currents in 
germanium and has been treated theoretically. 6 4 Dacey * has also in- 
vestigated the effect of the square root law upon the D(t) curves for the 
p-n-p structure of Section 4 and reports that the effects are so unfavor- 
able that no negative resistance is to be expected. 

The staueffekt opens the attractive possibility of making negative 
resistance devices in which the current decreases with increased dc 
voltage so that negative resistance will be exhibited over a wide fre- 
quency range. Unfortunately, when the boundary conditions are taken 




Fig. 5.1 — Qualitative representation of drift velocity versus field as affected 
by "staueffekt." 
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into account, it is found that a device in which most of the current flow 
occurs in a negative /z* region does not necessarily show a dc negative 
resistance characteristic. On the other hand, such a structure may have 
a very favorable I)(t) characteristic. 

We shall illustrate these conclusions by considering a (p + )p(p + ) 
structure having heavily doped ends, so that ohmic non-injecting con- 
tacts may be made to the ends. Fig. 5.2 shows the potential distribution 
and hole distribution for two cases of applied potential. The first case, 
represented in (a) and (b), corresponds to moderate fields such that the 
peak velocity u m is not reached. 

In the second case the voltage is so high that the average value of 
E = V/L exceeds the critical field E m . Under these conditions u is 
approximately equal to u m over a large part of the P-region and a sub- 
stantial portion of the voltage drop occurs near one end. An increase of 
the applied voltage occurs chiefly at this end with a small increase in 
current. Thus no negative dc resistance occurs. 

The above conclusions are reached by considering the differential 
equation for the space charge again as in Section 4 neglecting diffusion 
and starting with E = at the left edge of the P-region. This leads to 



dx = KdE/[(J/ii) - Ptt \, 
where we have introduced 



(5.3) 



Pa = -p f = -qN a (5.4) 

for the charge density of the acceptors. Evidently if J exceeds J m where 

Jm = p a u m , (5.5) 
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Fig. 5.2 — Hole density and potential distribution including influence of 
'staueffekt." 
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the denominator is everywhere positive and x increases monotonically 
with E. If ./ is only a little larger than J m , there will be a large region of 
x in which E is nearly equal to E m . Outside of this region, E increases 
much more rapidly with x. The relative scale of these distances may be 
estimated as follows: Suppose J is only a little larger than J m , and con- 
sider the situation where E is about twice E m so that u is about one half 
of u m . Then 



dx/dE = K/ Pa 



(5.6) 



Under these conditions an increase of E by an additional E m will require 
a distance 



Ax = KEJ Pa . 



(5.7) 



If this value is much smaller than L, then the situation represented in 
Fig. 5.2(c) and (d) will occur. On the other hand if L is smaller than 
Ax, the region of space charge and high field will extend throughout 
most of the structure. 

In any event equation (5.4) leads to positive resistance. This can be 
seen from the fact that increasing J always means a decrease in x for 
the same value of E and hence an increase in E at all values of x and 
thus an increase in voltage at any fixed value of x. 

The above conclusion that a positive dc resistance will be exhibited by 
a structure like that discussed above may also be reached by considering 
the transient response. The theory of Section 4 may be once at be applied 
to this case by simply taking account of the fact that v is negative for 
part of the structure and thus that 8E m increases with increasing s. 

In Fig. 5.3 we illustrate a structure to which these considerations may 
be relatively simply applied, at least in a limiting case. It consists of 
four layers, the two outer being p + as before. Space-charge limited emis- 
sion then enters the intrinsic layer which is of such a width that at its 
right hand boundary the electric field has a value E 3 that exceeds E m . 
At this point the hole space charge is 



pi = J/u z 



(5.8) 
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Fig. 5.3 — A structure having a region of uniform negative differential con- 
ductance. 
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where u 3 is u(E 3 ). In the P-region this space charge is compensated by- 
acceptors to produce a region of uniform field in which n* is negative. 
If the P-region is wide compared to the /-region, then the transit time 
through it will also be relatively large. As a consequence $Q will be 
transferred quickly into the P-region. From that time on 8E m > curves, like 
those of Fig. 4.3, will show an exponential increase with time and also 
with distance since for this case of constant u in the P-region, time and 
distance are linearly related. This will lead to a D{t) of the form 



D(t) = (u 3 /K)(S - t) exp | M * pz/K \ t, 



(5.9) 



where the absolute value signs emphasize that for this case of negative 
n* there is a build-up in time. This form of D is always convex upwards 
and, in fact, if 



S | M * pi/K | > 1, 



(5.10) 



it starts with a positive slope at s = so that the transient voltage 
actually builds up initially with time. 

Even an initially growing D(t) does not give a negative resistance at 
low frequencies, however. As shown in Section 2, the dc resistance is 
simply the integral under the D(t) curve and thus will still have a posi- 
tive value. 

5.2. Convergent Geometry 

It is possible to obtain marked improvement of the D(i) curves without 
the aid of the negative values of n*. This possibility is based upon con- 
vergent geometry. A possible case is illustrated in Fig. 5.4. In this case 
it is supposed that the field in the inner P-region is so large that a sub- 
stantial reduction in n* has occurred. As a consequence, the decay of field 
in this region is relatively slow. Furthermore, since both the dc and 




Fig. 5.4 ■ — A convergent flow structure. 



• 



826 THE BELL SYSTEM TECHNICAL JOURNAL, JULY 1954 

transient fields are high in this region, essentially because of the inverse 
square law, the principal contribution to D(t) comes from this region. 
These two factors — - relatively slow dielectric relaxation near the center 
and principal contribution to D(t) from near the center — combine to 
give a D(t) characteristic which holds up well until the pulse of injected 
holes reaches the inner region. This may result in a favorable convex 
upwards D(i) characteristic. 

ACKNOWLEDGEMENTS 

The writer is indebted to a number of his colleagues for helpful dis- 
cussions and to W. van Roosbroeck for the calculations for Fig. 4.4, and 
to R. C Prim for Fig. 4.5. 

REFERENCES 

1.1. W. Shockley and W. P. Mason, J. of Appl. Phys., 26, No. 5, p. 677, 1954 

3.1. W. Shockley, Electrons and Holes in Semiconductors, D. van Nostrand, i\ . Y., 

1950, p. 312. 

3.2. See Reference 3.1. 

3.3. W. Shockley, B. S. T. J., 28, p. 435, 1949, Section 2.4. 

3.4. See References 3.1 or 3.3. 

4 1 W. Shockley and R. C. Prim, Phys. Rev., 90, pp. 753-758, 1953. See also G. C. 
Dacey, Phys. Rev., 90, pp. 759-763, 1953, and W. Shockley, Proc. I.R.E., 
40, pp. 1289-1314, 1952. 

4.2. W. Shockley and R. C. Prim, Reference 4.1. 

4.3. E. J. Ryder, Phys. Rev., 90, pp. 766-769, 1953, references. 

4 4. See Shockley and Prim, Reference 4.1. 

5.1. This theory of mobility in higher fields has been published by H. Kromer, 

Zeits f. Physik 134, pp. 435-450, 1953. Kromer considers the theory in con- 
nection with the values of a in point contact transistors but does not ex- 
plore it as a power source. The present writer derived the same result in a 
more primitive form in 1948 as a potential means of obtaining high fre- 
quency power. However, experimental results by E. J. Ryder did not show 
evidence of this effect. The effect may possibly have occurred without being 
recognized because of the absence of an adequate appreciation of the im- 
portance of the boundary conditions discussed in this section. 

5.2. See Reference 4.3. 

5 3. G. C. Dacey, Phys. Rev., 90, pp. 759-763, 1953. 

5.4. W. Shockley, B. S. T. J., 30, pp. 990-1043, 1951. 

5.5. Personal communication. 



